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Abstract. We consider a two-phase problem for two incompressible, viscous 
and immiscible fluids which are separated by a sharp interface. The problem 
arises as a sharp interface limit of a diffuse interface model. We present results 
on local existence of strong solutions and on the long-time behavior of solutions 
which start close to an equilibrium. To be precise, we show that as time tends 
to infinity, the velocity field converges to zero and the interface converges to a 
sphere at an exponential rate. 
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1. Introduction 

We study the flow of two incompressible, viscous and immiscible fluids inside 
a bounded domain ft C R", n = 2,3. The fluids fill domains fl + (t) and Fl~(t), 
t > 0, respectively, with a common interface T(t) between both fluids. The flow is 
described in terms of the velocity v : (0, oo) x ft — > 1" and the pressure p: (0, oo) x 
D, — »■ M in both fluids in Eulerian coordinates. We assume the fluids to be of 
Newtonian type, i.e., the stress tensors of the fluids are of the form T(v,p) = 
2/u, ± Dv — pi in r2 ± (t) with constant viscosities fi^ > and 2Dv = \7v + Vv T . 
Moreover, we consider the case with surface tension at the interface. In this model 
the densities of the fluids are assumed to be the same and for simplicity set to 
one. For the evolution of the phases we take diffusional effects into account and 
consider a contribution to the flux that is proportional to the negative gradient of 
the chemical potential /i. Precise assumptions are made below. This is motivated 
e.g. from studies of spinodal decomposition in certain polymer mixtures, cf. |28j . 

To formulate our model we introduce some notation first. Denote by &T(t) the 
unit normal of T(t) that points outside fl + (t) and by V and H the normal velocity 
and scalar mean curvature of T(t) with respect to vv(t)- By [■] we denote the jump 
of a quantity across the interface in direction of Vr(t)i i- e -i 

x) = lim(/(a; + hv T n\) - f(x ~ hv T t t) )) for x e T(t). 

h^Q 
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Then our model is described by the following system 



dtv + v ■ Vu — div T(v,p) 


= 


in fi ± (t) for t > 0, 


(1.1) 


div v 


= 


in n ± (i) for t > 0, 


(1.2) 


mA/i 


= 


in fi ± (t) for t > 0, 


(1.3) 


~vr(t) ■ {T{v,p)\ 


= aHv r(t) 


on T(t) for t > 0, 


(1.4) 


V - u T (t) ■ v\ F (t) 


= -m\v T(t) 


• V/x] on r(t) for t > 0, 


(1.5) 


Mlr(t) 


= aH 


on r(t) for 4 > 0, 


(1.6) 



together with the initial and boundary conditions 



v\ au = on dn for t > 0, (1.7) 

uq, ■ mVfJ>\dn = on dfl for t > 0, (1.8) 

« + (0) = m, (1.9) 

u|t=o = hi ft, (1-10) 



where are gi ven initial data satisfying dfl^ fl c?r2 = and where er, m > are 

a surface tension and a mobility constant, respectively. Here and in the following 
it is assumed that v and [i do not jump across T(t), i.e., 



H = M =0 on r(t) for f > 0. 



Equations (|l.ip - (|1.2j) describe the conservation of linear momentum and mass in 
both fluids and (|1.4p is the balance of forces at the boundary. The equations for 
v are complemented by the non-slip condition (|1.7p at the boundary of CI. The 
conditions (|1.3|) . (|1.8|l describe together with ()1.5[) a continuity equation for the 
masses of the phases, and (|1.6|) relates the chemical potential fi to the I/2-gradient 
of the surface area, which is given by the mean curvature of the interface. 

For m = the velocity field v is independent of fx. In this case, (II. 5[) describes 
the usual kinematic condition that the interface is transported by the flow of the 
surrounding fluids and (|l.ip - (|1.10p reduces to the classical model of a two-phase 
Navier-Stokes flow as for example studied by Denisova and Solonnikov [10] and 
Kohne et al. [53], where short time existence of strong solutions is shown. On the 
other hand, if m > 0, the equations (|1.3p . (|1.6p . (|1.8[) with w = define the Mullins- 
Sekerka flow of a family of interfaces. This evolution describes the gradient flow for 
the surface area functional with respect to the i? _1 (f2) inner product. Therefore 
we will also call f| 1 . 1 [> - fjl . 10[) the Navier-Stokes/Mullins-Sekerka system. 

The motivation to consider (jl.ip - dl.lQI) with m > is twofold: First of all, 
the modified system gives a regularization of the classical model m = since the 
transport equation for the evolution of the interface is replaced by a third order 
parabolic evolution equation (cf. also the effect of m > in (|1.13p below). Secondly 
(|l.ip - (|1.10p appears as sharp interface limit of the following diffuse interface model, 
introduced by Hohenberg and Halperin [20] and rigorously derived by Gurtin et 
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al. pjj]: 

d t u + u-Vu-div(2^(c)Dw) + Vp = -£div(Vc(g>Vc) inf2x(0,oo), (1.11) 

divu = in ft x (0,oo), (1.12) 

3 f c + D-Vc = mA/j in x (0, oo), (1.13) 

H = e- 1 f'(c)-eAc inflx(0,oo), (1.14) 

v\ 9n = onffix (0,oo), (1.15) 

d n c\dn = d n /j,\du =0 on <9fi x (0, oo), (1.16) 

(u,c)| t=0 = (vo,co) in Q. (1-17) 

Here c is the concentration of one of the fluids, where we note that a partial mixing 
of both fluids is assumed in the model, and / is a suitable "double-well potential" 
e.g. f(c) = c 2 (l — c) 2 . Moreover, e > is a small parameter related to the interface 
thickness, \i is the so-called chemical potential and m > is the mobility. We 
refer to [2 [8] for some analytic results for this model and to [22l [18] for results 
for a non-Newtonian variant of this model. For some results on the sharp interface 
limit of (|l.llj) - (|1.17j) we refer to A. and Roger [5, Appendix] and A., Garcke, and 
Grim 0]. 

The purpose of this paper is to prove existence of strong solutions of (|1.1[) - (|1.10[) 
locally in time. Moreover, we will prove stability of spheres, which are equilibria 
for the systems. (More precisely, we show dynamic stability of the solutions v = 0, 
p,p = const., and J7 + (i) = Br(x) C for all t > 0.) Existence of weak solutions 
for large times and general initial data was shown in [5] . 

In the following we will assume that ft C R n , n = 2, 3, is a bounded domain with 
C 4 -boundary and that fj, ,m,a > are constants. One essential feature of (|l.ip ~ 
(|1.10|) is the coupling of lower order between the velocity field v and the chemical 
potential fi in equation (|1 .5() . Indeed, we will obtain functions in the regularity 



classes n 6 L P (J; W^{Q\T{-))) and 

v G Hl^L^Y) nL 2 (J;H 2 2 (n\T(-)T). 

Taking the trace to T(t) yields V/z|r G L p { J; Wp 1_1 ^ p (r(-))") and by complex in- 
terpolation and Sobolev embeddings we obtain 

v e Hi(J;L 2 (iir)nL 2 (J;Hl(Q\T(-)r) L q (J; W^(Sl\T{-)) n ), 

where q > p and p < 2(n + 2)/n. This shows that the trace 

v\ r eL q (J;W^(T(-)r) 

possesses more regularity with respect to time compared to V/x|r- We make essen- 
tial use of this fact by applying the following strategy for the proof of local-in-time 
wcll-posedness. After parameterizing the free interface T(t) via the Hanzawa trans- 
form by a height function h, the basic idea is to reduce (ll.lj) - (|1.10j) to a single 
equation for h. To this end we first assume that the interface, hence h, is given. 
Then we solve the (transformed) two-phase Navier-Stokes equations to obtain a 
solution operator v = SnsQi). Doing the same for the (transformed) two-phase 
Mullins-Sekerka equations, this yields a solution operator // = SMs(h)- Finally, 
we consider the transformed evolution equation (|1.5jl for the height function h and 
replace v and p, by SWs(/i) and Sms (h), respectively, to obtain a single equation for 
h. This quasilinear parabolic equation in turn can be solved by parabolic theory. 
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The only point one has to take care of is that the solution operator Sns in nonlocal 
in time and space. Therefore one has to deal with a parabolic equation with local 
leading part and lower order perturbations which are nonlocal (in time and space) . 
Having solved the single equation for h readily computes the velocity, the pressure 
and the chemical potential by the solution operators obtained before. 

Let us comment on the choice of an L 2 -setting for the Navier-Stokes part, while 
the equations for the height function h and the chemical potential \i are treated 
by an L p -theory, p > 2. One advantage is that the optimal regularity result for 
the two- phase Navier-Stokes equations with a given interface (see Theorem [AT]) is 
more or less easy to prove since it relies solely on resolvent estimates in L 2 . Another 
benefit is the reduction of the regularity of the initial velocity and the compatibility 
conditions at t = 0. For instance, if p = 2, then there is no compatibility condition 
for the initial value vo coming from the jump of the stress tensor, that is equation 

The structure of the paper is as follows: First we introduce some basic notation 
and auxiliary results in Section [2] Then we will prove that for a given sufficiently 
smooth interface T(t) the Navier-Stokes part of the system, i.e., (jl.lj) - (|1.2j) . (|1.4|) . 
(|1.7j) . (|1.10[) possesses for sufficiently small times a unique strong solution v in L 2 - 
Sobolev spaces, which are second order in space and first order in time. This result 
is proved using a coordinate transformation to the initial domains fij which goes 
back to Hanzawa and applying the contraction mapping principle. A key tool in 
our analysis will be a maximal L 2 -regularity result for the linearized Stokes system, 
which is proved in the appendix. Afterwards in Section [4] we prove that the full 
system possesses a strong solution locally in time for sufficiently smooth initial data 
by reducing the whole system to a single equation for the height function h (see 
above) . Then in Section [5] we prove stability of the stationary solutions that are 
given by v = 0, fi,p = const, and T(t) = dB r (xo) C fl and we show that (v(t), T(t)) 
converges to an equilibrium as t — >• oo at an exponential rate. 



2. Preliminaries 

2.1. Notation and Function Spaces. If A is a Banach space, r > 0, x e A, 
then Bx{x, r) denotes the (open) ball in A around x with radius r. We will often 
write simply B(x, r) instead of Bx(x, r) if X is well known from the context. 

The usual L p -Sobolev spaces are denoted by W k (Q) for k 6 No, 1 < p < oo, and 
H k (Q) = W$(n). Moreover W^ (n) and Hfi(Q) denote the closure of Cg°(f2) in 
Wp(£l), H k (fl), respectively. The vector-valued variants are denoted by Wp (f2; X) 
and H k (il; X), where A is a Banach space. The usual Besov spaces are denoted by 
B* iq (R n ), s G R, 1 < p,q < cx), cf. e.g. 0[36]. If ft C E n is a domain, B^ q {Cl) is 
defined by restriction of the elements of B^ g (K ra ) to ft, equipped with the quotient 
norm. We refer to [71 [33] for the standard results on interpolation of Besov spaces 
and Sobolev embeddings. We only note that (ft) and W k (O) are retracts of 
Bp q (W n ) and Wp(R n ), respectively, because of the extension operator constructed 
in Stein |35[ Chapter VI, Section 3.2] for bounded Lipschitz domains. In particular, 
we have 

«(0), < +1 (ft)) e , P = B k p + e (n) if ~ = i-^ + k e No, (2.1) 
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for all 6 G (0,1), cf. [361 Section 2.4.2 Theorem 1]. We also denote W k+8 {Q) = 
Bp+ e (fl) for fc £ No, 9 E (0, 1), 1 < p < oo. Furthermore, we define 

£(„)(«) = |/€L 2 (n) :^/(x)dx = o|, 

^(O) = {/eC -(fi)»:div/ = 0} i2(O) ". 

In order to derive some suitable estimates we will use vector- valued Besov spaces 
B s q ^(1; X), where s G (0, 1), 1 < q < oo, / is an interval, and X is a Banach space. 
They are defined as 

B q>oc (I;X) = {/ G L«(I;X) : \\f\\ B ^ { i;X) < oo} , 

\\f\\B-(I;X) = \\f\\L«(I;X)+ SUp \\A h f(t)\\ Lq{Ih;X ), 

0<h<l 

where A h f(t) = f(t + h) - /(f) and I h = {t E I : t + h E I}. Moreover, we set 
C S (I;X) = S^ i00 (J;X), s E (0, 1). Now let X ,X t be two Banach spaces. Using 

/(*) — /( s ) = Is lif( T ) di~ it is easy to show that for 1 < q < qi < oo 

W^(I;X 1 )nL*>{I;X )^B e q>oo (I;X e ), - = — + -, (2.2) 

q qo qi 

where 6 G (0, 1) and Xg = (X ,Xi)[ e ] or X = (X ,Xi)e )r , 1 < r < oo. Further- 
more, 

B l,oc{^ x ) ^C e -i(I;X) foraflO<0<l,l<g<oo with - - > 0, (2.3) 

cf. e.g. [35]. Furthermore, for s G (0,1) we define H s (0,T;X) = B$ 2 (0,T;X), 
where / E 5f j2 (0, T; X) if and only if / G L 2 (0, T; X) and 

rT ' T ||/W-/(r)» a 



\B° t2 (0,T;X) ~ ll/lll='(0,T;X) ^ / q J q \ f _ r | 



l -^dtdT < OO. 



In the following we will use that 



./0 

for all < s < s' < 1, which implies 

II/II^(0,t ; x) < C s , s ,T?\\f\\ c ,, {%T] . x) for all / G C s '([0,T];X) (2.4) 

provided that < s < s' < 1, < T < 1. 

Furthermore, we note that the space of bounded fc-times continuously diffcr- 
entiable functions /: U C X — >• Y with bounded derivatives are denoted by 
BC k (U;Y), where X, Y are Banach spaces and U is an open set. Moreover, 
f E C k (U;Y) if for every x E U there is some neighborhood V of x such that 
f\ v EBC k (V;Y). 

We will frequently use the following multiplication result for Besov spaces: 

ll/flK lM(niM) <c;..*«ll/K llfl llfllln5 >fl ( 2 -5) 
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for all / G Bp i qi (M n ),g G ^ g2 (R") provided that 1 < p < p± < oo, 1 < q h q 2 < 
00 ' r > p7' and 

-r + n(i + i-l) < s < r, 

\pi p / + — ' 

cf. [23 Theorem 6.6]. Since W°{R n ) = B s pp {R n ) for every s G (0,oo) \ N, this 
implies that 

||/3lk-(R») < C s , p \\f\\ Wsm \\g\\ Ws(Mn) for all /,g G % S (K") (2.6) 

provided that s — ^>0, l<p< oo. Concerning composition operators, we note 
that 

G(f) e ^ 9 (]R n ) for all G G C°°(R) with G(0) = 0, / G B^ g (M n ) (2.7) 

provided that again s — ^ > 0, 1 < p, q < oo. This implies that J -1 G B pq (Q) 
for all / G J3p j9 (f2) such that |/| > c > if is a bounded Lipschitz domain. 
Moreover, the mappings / h> G(/) is bounded on B pq (W l ) under the previous 
conditions. We refer to Runst |29j for an overview, further results, and references. 
Furthermore, using the boundedness of / i— > G(/) one can easily derive that 

G(.)eC l (B p s i? (l")) 

for any G G G°°(R) with G(0) = 0. To this end one uses 

G(f(x) + h{x)) = G{f(x)) + G'{f(x)) + f G"(f(x) + th{x)) dt h{x) 2 

Jo 

together with ([22]) and the fact that (G"(/ + th)) te[0A] is bounded in B» q (R n ). 

Finally, by standard methods these results directly carry over to W p (T,), B pq (T,) 
if E is an n-dimensional smooth compact manifold. Then G(0) = is no longer 
required since constant functions are in B p 9 (E). 

2.2. Coordinate Transformation and Linearized Curvature Operator. 

In the following let E C SI be a smooth, oriented, compact and (n — l)-dimensional 
(reference) manifold with normal vector field i>s. Moreover, for a given measurable 
"height function" h : £ -> K let 

6 h : £->M n : x i-> x + h(x)u s (x). 

Then 9h is injective provided that ||/i||loo < a for some sufficiently small a > 0, 
where a depends on the maximal curvature of E. Moreover, we choose a so small 
that 3a < dist(£, <9f2). Then the so-called Hanzawa transformation is defined as 

G h (x,t) = x + X (d^(x)/4a)h(t,U(x))^(U(x)), (2.8) 

where ds is the signed distance function with respect to E, II(a;) is the orthogonal 
projection onto E, \ G G°°(R) such that x(s) = 1 for \s\ < | and x(s) = for 
\s\ > | as well as |x'( s )l < 4 for all s G M, and ||ft||L°° < cl- It is well-known 
that 9/j(.,t): f2 — >• fi is a G 1 -diffeomorphism. Hence F^ := 0^(E) = 6>/j(E) is an 
oriented, compact G fe -manifold if h G G fc (E) with ||/j||z,°o(jq < a- 
For the following let 

U = j/iGWp 4 ~*(£) : \\h\\ L oo <a|, (2.9) 
Ei,t = i p (0,T;Wp"^(E))n W£(0,T; Wp 1_ *(E)) 
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where 3 < p < 2 , < T < oo. Furthermore, let 

K(h):=H h o9 h , (2.10) 

where Hh : r n — > R denotes the mean curvature of = 0/j(E), i.e., it is the sum 
of all principal curvatures. 

Lemma 2.1. Let 3 < p < ^"^ ' and U C W p p (E) be as above. Then there are 
functions 

P G C ,1 (W,£(W^"'(S),W^"^(S))), Q G C^W, WJ"»(E)) 

swc/i i/iai 

K{p) = P(p)p + Q(p) for allpeUD Wp 4 ~ p (E). 
Moreover, ifE, — Sr '■= dBn(0), then 

1 /n-1 \ 
DiT(O) = L> := 1?^ := -— -1—+A Sr ). (2.11) 

Proof. The proof follows essentially from the proof of [HI Lemma 3.1] and [121 
Remark 3.2 a.]. To this end let {(Ui,ipi) : 1 < I < L} be a localization system for 
E, i.e., E = Uj_j. Ui and <yS; : (— a,a)™ _1 —> Ui is a smooth local parametrization of 
C/j for alH = 1, . . . , L. Moreover, let s — (si, . . . , s„_i) be the local coordinates of 
Ui with respect to this parametrization and 

pi(s) :=p{(pi{a)), X t (a,r) := X(ipi(a),r), (s,r) G (-a,a) n 

be the local representations of p, X, where I: Ex (—a, a) — >• R" with X(s,r) = 

s + r^ s (s) and p € U C W p P (E). Then it follows from [HI Equations (3.4), 
(3.5), Remark 3.2 a.] that K(p) — P(p)p + Q(p), where P(p),Q(p) have the local 
representations 

1 / n_1 \ 1 

p i{p) = — j E Pok(p)ds 3 d Sk +22P i (P) ds * ' Q'O 3 ) = 
71 y,*=i »=i / 

where 

1 / " -1 

p*o>) = 73 - l > Jk (p)+ E ^'^(fR^ 



1 / n— l n— l n— l 

p*0>) = js U E ^" r lfe + E ^ ki r%d slP + Y 2w km r nk d SmP 

P y ,j',fe=l j,i=l fe,m=l 

n-1 \ 

- 7,fc,Z,m=l / 



P 3,fc=l 



\ 



1 + X w3 k {p)d S] pd Sk p, 



n-1 



T )k (P) = E wm (<°K X-d Sm X\ Ms)) , i^n, 

r "fe(p) = d S] d Sh X ■ d Sn X\ (sAs)) , w jk (p)(s) = d S] X ■ d Sk X\ (sAs)) 
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and (w j7£ (p)(s))" fe i 1 is the inverse of (w jk (p)(s))l k i 1 . 

4-A 

Since E is smooth, X and d Sj X -d 8 .X are smooth. Therefore Wjk(p) £ W p p (E) 
because of (|2.7p . Since det((w J fe)™^ 1 ) > c > by construction, we obtain 

w ]k (p) £ W^"* (E) for all = 1, . . . , n - 1 because of (J2JJ). 
Moreover, 9 Sj /5 G W p p (E) and therefore 

"- 1 4 
£ w^ k (p)d Sj pd Sk p £ Wp~^{Y<) 

due to (|2.6p . Using p. 71) again, we obtain Z p £ W p p (E). Proceeding this way, we 

finally obtain that Pjk(p),Pi(p),q(p) £ W p " (E) for all p £U. Now (|2.5[) implies 
that 

2-i < Cp||a|| 3 _4 ||u|| 2 _i 
w p p (s) F w p "(E) w p 

for all a £ Wp~ 5 (E),u G Wp~*(£). Hence 

P € C 1 (W, £(W P ~ ' (E), W p 2 ~» (E)), 

since the operators are compositions of C 1 -mappings. Moreover, (|2.11[) follows 
directly from the observations in the proof of |12[ Lemma 3.1]. □ 

Corollary 2.2. Let K be as in ||2J0]) . Then 

K £ C^Elt n U- Hi (0, T; Z, 2 (E)) n L 2 (0, T; #3 (£))). 

Moreover, for every e > 0, < To < oo i/iere is some C > suc/i i/iai 

^^BC 1 (E lir nW e ;Hi(0,T;L 2 (E))nL2(0 : T;Hi(E))) — ^ 

for all < T < To, where U e = {a £ U : ||a||i«>(E) < a — e}. 
Proof. We use that 

X(/i)= ^ a a (x,h,\7 s h)dfh 

\a\<2 

for all ft G C 2 (E), where a Q : E x R x M" _1 -> K is smooth. Since 

Ei, T ^ bIooCO.T;^ "*(£)) nS| >oo (0,T; Wp 3 "^(E)) 
due to (|2T2|) and 

4,oo(0,T;W^~-(E))^C7*([0,r];C (S)) 
due to (I2.3[) and p > 3, we conclude that 

o a (a;,ft, V 5 /i) G C*([0,T];C Q (£)) for all h £ E^ T nU 
and for all \a\ < 2. Moreover, the mapping 

WnEi, T 9 h^f a a (x,h,V s h) £ C^([0,T]; C*°(E)) 
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is C 1 since a a are smooth. Furthermore, we conclude that 

<ai| aQ (,,^v s ^>|| sL(oT;ip(s)) 

<C^a a (x,h,V s h)\\ i \\v\\ , t _i 

< all«a(^^V^)|| a([oT];CO(s)) || V || EliT 

for all \a\ < 2, v £ E^t, ft £ E^t nW e , £ > 0. Since multiplication is smooth (if 
bounded), it follows that 

K G BC^Ei.t nW £ ;ff' (0, T; L 2 (E))) 

for any e > 0. Finally, we use that a a (x, ft, V s h) G BUC([0, T]; C 1 (S)) and 

E liT nW e 9 ft h> o a (£c,ft, V s ft) e BZ/CftQjT];^ 1 ^)) 

is in C 1 with bounded derivative. Hence 

a a (x, ft, V s ft)V> G L p (0, T; Wp - * (E)) i 2 (0, T; ffs (£)) 

for every h E U £ tl E^t, e > and the mapping ft i— > -fT(ft) is in BC 1 with respect 
to the corresponding spaces. Altogether we have proved the corollary. □ 

3. Two-Phase Navier-Stokes System for given Interface 

In this section we assume that the family of interfaces {T(t)} t >o is known and 
we will solve the system (fTTjl . (jL~2]) . (fL~4|) . (fTT)) . (fTTTU)l together with the jump 
condition [u] =0. 

For the following let E C be a smooth compact (n — l)-dimensional reference 
manifold as in the previous section. Moreover, we assume that there is a domain 
fij CC f2 such that E = 8Qq . Moreover, we assume that 

r(t) = {x + h(t, x)isx(x) : x G E} =: r ft(t) 

for some ft eW flEu, where 

Ei, T :=% 1 (J;X )nL p (J,X 1 ), 

J = [0,T], and 

= Wp~'(Ti), Xi = VFp"^(E) 

for p > max( :s ±^,3) = 3, n = 2,3, and ^s(^) is the exterior normal on dti^ = E. 
Here W is as in ([279]) . 

For given ft G E^t let ft = -Eft G E^t, where 

£: Ei, T -> 1i,t := H^(J; W£ (!!„)) n L p (J; VF p 4 (E Q )) 

is a continuous extension operator and E a = {a; G fi : dist(ic, E) < a}. Then by 
Lion's trace method of real interpolation, we have 

E hT ^BUC{[0,T};X, ( ), X, ( = Wp~' (£«) ^ C 2 (ET) (3.1) 
since p > 2i j^- Moreover, if we equip and Ex t with the norms 

IMk T = ||u|| i-i 4-i +IKo)||x 7 , 

p (E))nip(J,Wi, P ( s )) 

IMIl 1T = ll"lk p i(,/;H' p 1 (s a ))nL p (x^(s a )) + IK0)|| x , 



10 



HELMUT ABELS AND MATHIAS WILKE 



then the operator norm of the embedding ()3.1[) is bounded in T > 0. Additionally, 
we have 

Ei, T ->C 1 -»([0,r|;W r p 1 (E o )). 
Interpolation with (|3.ip implies 



E 1)T C T (\0,T]-B p y^ a )) C r ([0,T];C 2 (E a )) 

for some r > since p > ^jF*. Here again all operator norms of the embeddings 
are bounded in T > 0. We will need the following technical lemma: 

Lemma 3.1. For every e > the extension operator E above can be chosen such 
that for every < T < oo 

sup \\h(t,U(-)) ~ Eh(t r )\\ c i^ a) < e\\h\\ El T . 

0<t<T 

Proof. First of all, since Eh(t, x) = h(t, x) for all x G E, t G [0, T], 

sup ||/i(i,n(-))-^/i(*,-)llcHs a ,) <«' sup \\Eh(tr)\\c*p u ,)<Ca'\\h\\^ T 

0<t<T Q<t<T 

for any < a' < a, where C is independent of < T < oo. Hence, if, for given 
e > 0, a' is chosen sufficiently small, we have 

sup \\h(t,U(-))~Eh(t,-)\\ c i^ al) <e\\h\\ ElT . (3.2) 

0<t<T 

If we now define E' : E^t — » E^t by 

(E'h)(t,x) = (Eh) (t,H(x) + ^-d s (x)i/ E (n(x))^ for all x G E a , t G [0,T], 

then E' : E^t — > ^i,t is an extension operator, which satisfies the statement of the 
lemma. □ 

For technical reasons, we modify the Hansawa transformation 0^ to 
@h(x,t) = x + x(di:(x)/a)h(t,x)i^(Il(x)), 
where h = Eh G E\ t T is the extension of h to £1 as above. Then 

\\e h (.,t)-e h (.,t)\\ cim < c\\h(.,t)-h(n(.) 7 t)\\cHx a) 

for all < t < T, where C is independent of h and < T < oo. If we now choose 
e > in p. 21) sufficiently small, 0/j(.,t): fl — > is again a C 1 -diffeomorphism for 
every < t < T . This can be shown by applying the contraction mapping principle 
to 



,T = 9- 1 (e h (x) - Q h (x) + y) 



for given y G f2, which is equivalent to ®h{ x ) — V- Moreover, ©^(E, t) = 0fc(E, t) 
T(t) for all < t < T. 
Now let 

F M = e h (.,t)oe h (. ) o)- 1 . 
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Then F h ,t : -> with F/, )t (n^) = fi^i) and F h ,t(T ) = T(t), where T = T(0) = 

Ol+(0). Moreover, F h = (F h , t ) te[ o, T ] G BZ7C([0, T]; Wp~^ (Q,)) n W^(0, T; W^(O)) 
and 

H-F/jj - -F 1 / l2 || c , T([0 T j. c , 2(n)) < C||/ii - /i2||e 1iT ) ( 3 - 3 ) 

llFij - F/j 2 ||wi(0,T;W£(n)) < C||/li — /12||Ki jT ( 3 - 4 ) 

for all ||/ij||Ei t < R, j ' == lj 2, where C is independent of /ij and < T < oo. Since 
Fh : o = Ida for all h G Ei,t, (|3-3p implies 



H-F/i! - -P)s 2 ]| b£/C ([ 0)T ]. C 2(q)) < CT T ||/ii -/i2||Ei, T - (3.5) 
Now we consider 

d t v + v Vv- n ± Av + Vp = in 0^(4), t G (0,T), 

divv = in ± (t),t G (0,T), 

M=0 on r(t),f G (0, T), 

• T(v,p)] = aH r{t) v m on T(t),t G (0,T), 
u| 9 n = on an,te(0,T), 

w|t=o=u onO ± (t),iG(0,T). 

Defining 

= v(F t>h (x),t), q(x,t) = p(F t> h(x),t), 
the latter system can be transformed to 



- ( u ± Au + Vq 


= a ± (/i; D x )(u,q) + d t F h ■ V h u - u • V h u 


in Q T , 


(3.6) 


div u 


= Tr((J - A(/i))Vu) =: £f(ft)u 


in Q T , 


(3.7) 


["1 


= 


on r 0> T, 


(3.8) 


\v To ■ T(u, q)j 


= t(h;D x )(u,q)+aH h 


on r ,T, 


(3.9) 


u\on 


= 


on 


(3.10) 


U {=0 


= «o 


On fig, 


(3.11) 



where Q% = (0,T)xOJ, ^„ = O\(r UfiJ), r , T = (0,T)xr o , dO T = {0,T)xdil. 
Here 



a ± (ft;D a! )(tt,g) = /i div^(V^u) - ,(/ div Vu + (V - V/,)g, 

V h = A(ft)V, div h u - Tr(Vhw), - DFf^ , v h - 

\ A \ h ) v v a \ 

t(h,D x )(u,q) = [(vr a - v h ) ■ (2fi ± Du - ql) + 1v h • sym(Vw - V h u)j, 
H h {x) = H T{t) {F h ,t{x))v T{t) {F h , t {x)) for all x G T . 

In the following let Yr = Y£ X 7|, where 

Yj* = £ B[/C([0,r];iJ 1 (t])' 1 ) ntf^T;^^) : € L 2 (0,T;ff 2 (fi±)")} 
Y? = [q G L 2 (0,T;L 2 , (0) (n)) : V g | ± G L a ((0,T) x fi±)"} . 
The main result of this section is: 
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Theorem 3.2. Let R > 0, h G U. Then there is some T = T (R) > such that 
for every < T < To and h G E^tHW with h\t=o = ho and vq G Hq (f2) n nl2,rr(fi), 
n = 2,3, with max{||/i||£ 1 TJ ||^o||.ffi(n)} ^ R there is a unique solution (u,p) =: 
J"r(/i,«o) S It °/ (|3.6p - (|3.1ip . Moreover, for every e > 



JV e BCH^b x B H i (0, R); Y T ), 



where 



A £iR = {he B ElT (0,R):h(0) = h Q , sup ||/i(t)||ioc (s) 

0<«<T 



< a — £ 



We can formulate ()3.6|) - (|3.11[) as an abstract fixed-point equation 

Lw = G(w; h, vq) in Zt 

for w G ir, where 

/d t u~ n ± Au + Vq\ 
div u 

U|t=0 



(3.12) 



L{u,q) 



\ 



J 



G(u,q;h,v ) 



/a ± (h; D x )u + d t F h ■ X7 h u - u ■ Vhu\ 
g(h)u- j^ J n g(h)udx 

t(h;D x )(u,q)+aH h 
V vo J 



for all w = (u, q) G Y T , where Zt = Z T X Z T X Z T X Z, 

Z T = L 2 ((0,T) xQ ) n , 4 = ff l (n) n nL v (O), 

Z| = L 2 (0,T; J ff ( 1 0) (r!o))ni/ 1 (0,T;ff ( - ) 1 (r! )), 

z| = L 2 (o,T;fl 2 »(ror)nfri(o s T;L 2 (ro) n ) ) 

and Z| = H%(Q,) n n L 2 . a {VL). Here ff^fio) = ^(^o) H £ 2 ,(o)(^) is normed by 
||V.|| i2 ,ff ( - ) 1 (fi) = ( J ff ( 1 0) (fi))'. 

First of all, let us note that f|3 . 1 2[) implies (|3.6I) - (I3.11I) except that (|3.7[) is replaced 

by 

div u = g(h)u - -!- J g(h) dx. 

But the latter equation implies (|3.7p . which can be seen as follows: Let K(t) = 
1^1 J a g(h(x,t))dx and u(af,i) = u^ 1 (a, «),*)■ Then v(t) G £#(0) for all t G 
(0, T) and therefore 

divv(x,t)dx= / Tr(A(/i(a;,t))Vit(a;,t))detDi ;, / l (af ) t)eiE 



= K(t) I det DF h (x,t)dx 
Jn 

for all i G [0, T]. Since the last integral is positive, we obtain K(t) = for all 

te(o,T). 

Lemma 3.3. Le£ i? > 0, e > 0, and let Yt,Zt, ho be as above. Moreover, let 
A e ,R = \ h G Ei, T : sup ||/i(t)||z,-(s) < o-£,h(0) = ft , INkr < #f • 

L 0<t<T J 
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Then there is some To > such that for every < T < To the mapping G defined 
above is well-defined and 

G e &{B Yt $,R') x A £<R x B H i nL3 jO,R'); Z T ). 
Moreover, there are some C, a > such that 

\\G( Wl ;h,v Q ) -G(w 2 ;h,v )\\z T < CT a \\ Wl - w 2 \\ Yt 
for every u>i,w 2 E B Yt (0,R'), < T < T a , h e A S>R , and v G B H i nL . 2 a (Q,R). 

Proof. First of all, because of (|3.3[) , for any e > there are some C, T > such 
that 

|| Ar-FX - id Hbc/ccio.t]^ 1 ^)) - ^II^IIei.t 

for all < T < T , ||/i||ei t < R- Hence F t) /j : O — >• f2 is a C 2 -diffeomorphism and 
D x Fh is invertible with uniformly bounded inverse for these h, T. Since matrix 
inversion is smooth on the set of invertible matrices, 

A{h) = DF^ T G C T ([0,T];C\n)) 

for some r>0if||/i||E lr <-R- Moreover, interpolation of (|3.3p and (|3.4j) yields 

DF h G C^ + ^([0,T];C°(n)) 

due to W^(0,T;Z) ^ C 1 ^ ([0, T]; X), where the operator norm of the latter 
embedding is bounded in < T < oo if VFp (0, T; X) is normed by 

\\f\\w i{ O t T;X) := \\(f,f)\\L p{ 0,T;X) + ||/(0)||x. 

Here we have also used that ||/|| cl(S) < C\\f\\ J 0(n) ||/|| * a(n) and W£(fi) C°(fi). 
Hence 

= DF f ; T G C^* + *([0,T];C°(n)). 

Furthermore, 

A e SC^l^^O^RjjX) with (3.13) 

X = C r ([0, T]; C* 1 ^)) n W£(0, T; L p (Jl)) n C* _ & + * ([0, T]; O (fi)) (3.14) 

again since matrix inversion is smooth. 

Using the above observations, one easily obtains 

||(Vh - v)/|| L2(0T;iffc(a ± )) < cr" r ||/i||E 1(X ||/|| i3 ( 0)T . J jfe+i (n ±)) 

for all / 6 L 2 (0,T;H k (Q t )), k = 0,1, \\h\\ El T < R. From this estimate, one derives 

||a ± (/i;L> :r )(M,g)|| i2((0iT)><n ± ) < CT T \\h\\ El T \\(u, q)\\ YT , 

\\9( h h\\ L2( o,T;Hi(n±)) ^ CTT H' l llEi,Tll«ll£ a (o,r ! fl'(n±))' 
W v - DF h Vu \\ L2 ((o,T)xnt) ^ crr *H ,, ll2i»(o,T i wj(n±})ll u llL <x> (o 1 TiWj(n±)) 

< CTi\\v\\ n \\u\\ Y i 

\\d t F h ■ Vu|| ia((0)r)xn ±j = || (d t F h - d t F ) ■ VuHi^^rjxn*) 

< CT^p\\h\\ ElT \\u\\ Loo( o. T . H i {n)) 

< CT^\\h\\ EliT \\u\\ Y ^ 
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where we have used (13. 4|) for the last estimate. Moreover, 

\\(t(h,D x )u\\ Z 2 < CT a \\h\\ El . T \\(u,q)\\ YT (3.15) 

for some a > can be proved in the same way as in [TJ Proof of Lemma 4.3]. In 
order to estimate g(h)u £ -ff 1 (0, T; Hz} (f2)), we use that 

(g(h)u,<ph = -(u,div((J- A(hf)<p)) n for all ^ 6 Hf 0) (Sl). 
Therefore we obtain for all <p € H^(fl) with || V(^||i 2 (n) = 1 

^(ff(%^)n = -(^,div((/-A(^) T )^)) fi -(Vu,(5M(/i) T )^)fi 
= (F 1 (i) ) ^) + (F 2 (t),^), 

where 

|(a tU (i),div((/-A(^(i)) T )^))a| 

< c||a tM (t)iu 2(f , ) ||/-^W T || ioo(0 , T;C1( n ) ) 

< C*T r ||<9 t u(t)|| i2(n) ||/i|| CT([0T;C , 1(r?)) 

and 

\(\7u(t),(d t (I-A(h(t))) T ^) n \ 

< C\\u\\ L ^ (0 ,T-,H^Q))\\\\dtE(h(t))\\ w i {Q) 

< C\\u\\ n \\d t A(h(t))\\ w}m 
for all t G (0,T). Hence 

II F i|Il 2 (0,T;H ( -J) ^ C ' TT |l 5 t M ll£2(f2x(0,T))ll^llc-([0,T];C 1 (n)) 

11^11^(0,^) < C\\u\\ Y ,\\d t A{h{t))\W ,T,wl) <CTW|| u || y i||/i|| Elt:r . 
and therefore 

\\d t g(h)u\\ L2{0T . H - Q]) < C(R)T™^-^\\h\\ El ^u\\ Y r 

for all h G Ei.t with ||/i||ei t < R- Here we have used that A € _BC 1 (j4 £j ij; X), 
where X is as in (|3.14p . 

Finally, it remains to estimate the term H^. To this end we use that 

H h = (K(h)oO£) Vh 

where 9h ■= S/i(-,0)|e: E — > T bijectively. Here 

K G BC X {A Z<R ; Hi (0, T; L 2 (E)) n L 2 (0, T; ff* (E))) 

because of Corollary [221 Since 6>/j G C 2 (E) n is independent of t and /i, the same is 
true for K^oO^ 1 with E replaced by IV Because of p,13[) . we have for H(h) := H) x 
for all h G A e ^R 

H G BC^Be^R); Hi (0, T; L a (T )) n L 2 (0, T; JT* (T ))). 

Altogether, since all terms in G are linear or bilinear in (u,q) and A(h), these 
considerations imply that G G BC 1 (Sy T (0, i?) x A s> r x B ff i(s))n(0, i?); Z T ) and 

\\G(wi;h,v ) - G(w 2 ;h,v )\\z T < CT a '\\ Wl - w 2 \\y t (3.16) 
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for all Wj = (uj,qj) £ Yt with ||iOj||y T < R, h £ A S<R , vq € 5z 4 (0, R) and 
< T < T for some a' > 0. □ 

Proof of Theorem 13. 2t Let e > 0. Using Lemma 13.31 and choosing To > 
sufficiently small, 

L- X G{.; h, v ) : B Yt (0, R') B Yt (0, R>) 

becomes a contraction and is invertible if h 6 A e ^ and ||wollir 1 (n) ^ -R; where 

# = 2sup{||L- 1 G , (0;/i, Uo )||y T : h S A £)R) ||«o|| H1(n ) < fl} • 

Hence for every (/i, «o) S ^e,_R x Bz4.(0,R) there is a unique «; =: FriKvo) G 
B Yt (0,R') such that 

u> = L~ 1 G(w; h, vo). 

Moreover, (|3.16p implies 

||T-^G(i<;;Mo)IU(z r ) <CT"' < 1 

for all tOj = (itj.gj) e T T with ||w,-||y T < i?, (7i,«o) € x B H i nL2 a (0, i?), and 
< T < To if To is sufficiently small. Hence we can apply the implicit function 
theorem to 

F(w; h, vo) — w — L^ 1 G(w; h, vq) — 

and conclude that 

F T £ BC 1 (A e<R x BzA0,R);B Yt (.0,R')) 
since D w F(w; h, vo) is invertible for all w £ B Yt (0, R'), h £ A e ^ R ,vo £ B Z i(0,R). ■ 

Finally we obtain that the mapping h h4 {yh ■ u) o (0h|t=o)|E satisfies the condi- 
tions to apply the general result of [6]: 

Corollary 3.4. Let R,e > 0, T = T (R) > 0, A e , R , and Ft be as in Theorem [KR 
For every h £ A £ ^ R , v £ i?o(^)™ n L 2jCT (fi) with \\vq\\ h i^ < R let 

g T (h;v ) := {v h ■ u) o (8h|t=o)|E, 

where (u,p) = FT{h,vo), < T < To. Then there is some q > p such that 
Qt € C X (A^ R x B H i nL2 a (0,R);L q {0,T;X )). Moreover, if /n| [0>T '] = M[o,t'] 
/or some < T' < T, then GtQii; fo)|[o,r'l = GtQ^'i w o)|[o,T']> is., the mapping 
h i-> GT^h; vq) is a Volterra map in the sense of [6]. 

Proof. First let n = 3. Then by interpolation 

ffV; L 2 (0)) n T 2 (J; ff 2 (fi)) L 4 (J; F*(fi)) n T 2 (J; H 2 (ty) 
L q (J;H s (n)) L q {J-Wl{tt)) 

where 

since 3 < p < ^ and n — 3. If n = 2, we use that 

H\J- La(fi)) n T 2 (J; ff 2 (fi)) ^ L 4 (J; Wtf (fi)) -4 L 4 (J; W^(fi)). 
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Hence Tt e C 1 (A £ ,r x B H i(0,R)-L q (J; W£(to)) for some q > p. The rest of the 
first statement follows from the trace theorem, the fact that 

h^v h € C^Be, T (0, R), BUC^T^C 1 ^))), 

and that Qh\xe's,t=o' E — » To is a C 2 -diffeomorphism. 

Finally, the Volterra property follows easily from the fact that the solution of 
(nO)|) - ([gTTj) on a time interval (0,T) is also a solution of (l3ToV ([3"T[T) on (0,T') for 
any < T' < T (after restriction) and the uniqueness of the solution. □ 

4. Local Well-Posedness 

In this section we show that the system (ll.ll) - (|1.10l) admits a unique local-in-time 
solution by reducing the whole system to a single quasilinear evolution equation for 
the height function h. For this purpose we use the solution operator obtained in the 
previous section and the solution operator for the (transformed) chemical potential 
coming from (I4.6[l - (|4.8p . 

We transform (TOl) . ([i~5j) . (TLB) and (Ti~8l) to the fixed domain fl\E, with Sell 
as in the previous section, by means of the Hanzawa transform. This yields 

mA h r) = in n T \Y, T , (4.1) 

d t h- {vh ■ u) o (6fc| t =o)|s = -m\vh ■ V^r?] on E T , (4.2) 

T)\ s = aK(h) on E T , (4.3) 

m ■ mV?7 = on <9f2 T , (4.4) 

h\ t =o = h on E, (4.5) 

where Et := (0,T) x E, r](t,x) :— fj,(F t ,h(x),t) and K(h) denotes the transformed 
mean curvature operator. Assume that we already know a solution (u, h) £ x 
E^t. Then we may use Corollary 13.41 to write (y^ ■ u) o (8/j|t = o)|s = 5t(^;^o)- 
Consider the elliptic (time dependent) problem 

mA h r] = in ft T \E T , (4.6) 

ri\s = <rK(h) onE T , (4.7) 
uq, ■ mVr? = on dfl T - (4.8) 

If 

h e E^ T := {h e Ei, T : H[0,T]) e U}, 

where i7 C Wp 4_4/p (E) is a sufficiently small neighborhood of zero, then (|4.6|l - 
gU) admits a unique solution 77 =: S(h)K(h) € L p (0, T; Wp 2 (fi\E)). Defining 
B(h)r) :— m\vh ■ V^r/] we may reduce (|4.1j) - (|4.5D to a single equation for h: 

d t h + B(h)S{h)K(h) = g T (h: u ) on E T , h(0) = h on E. 

Employing the decomposition from Lemma 12.11 we may write 

d t h + A{h)h = F hT {h) + F 2 {h) on E T , h(0) = h on E, (4.9) 

where A(h) := B '(h) S '(h) P (h) , F hT (h) := g T {h;v ) and F 2 (h) := -B(h)S(h)Q(h). 
Note that A and F 2 are nonlocal in x but local in t, whereas F\ t T is nonlocal 
operator in t and x, but it has the Volterra property with respect to t. Firstly we 
show that F 2 £ C l (U; V^p^ 1/p (E)). By LemmaOwe have 

QeC 1 (I7;^- 1/, (E)). 
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Next we show that S G C^U; £{Wp~ 1/p (T,), W|(fi\£))). Writing 

n—1 n—1 

A h = J2 "ikOjdk + Y, a J d S 

j,fc=l j=l 

with coefficients 

a 'jk = a jk(x, h, Vh, V 2 h), aj = a,j(x, h, Vh, V 2 h), 
depending smoothly on (x, h, Vh, V 2 h), it is not hard to see that 

a%(-,h,Vh,V 2 h),a^(;h,Vh,V 2 h) G BUC(n\Y,) 2 , 
for all h £U. Here we used the embedding 

VK 2 ~ 4/p (S) C(£) 
whenever p > (n+ 3)/2. This in turn yields that 

h m> A/j 6 C 1 ([/; £(Wp (f2\E), 

We can now write 

S(% = (Ah,7,7JV,an) -1 (0,ff,0) 

for some function g G Wp 2 Here 7 denotes the trace operator to £ and jN.dn 

stands for the Neumann derivative on d£l. Since the mapping h H> (Aft, 7, 7^,90) 
belongs to 

and inversion is smooth, we may conclude that 

Finally, we show that B G G l {U; C{W 2 {VL\T), Wp~ 1/p (£))). We may write B(h) = 

Y^jZi bjjdj, where the coefficients bj — bj(x, h, Vh) depend smoothly on (x, h, Vh). 

This yields bj{-,h,Vh) G C l (T), for each h e U since Wp 3_4/p (£) C^E) for 
p > (n + 3)/2. It follows readily that 

h h+ 6^ G C^AW^E), W^-Vf (S))). 

Summarizing we have shown that 

F a GC 1 (l7;W p 1 - 1 /P(E)) J 

hence the desired assertion. 

Concerning the mapping /i 1— ^ A{h), we would like to show that 

h^A(h) g C l (U;C(W^(J:),W^(E)))- 
But this is an immediate consequence of Lemma 12.11 since 

h ^ P(h) G C\U; C{Xi, w 2 - 1 ^^))). 

It has been shown in [37J Proof of Theorem 4.1] that .4.(0) has the property of 
maximal L p -regularity in X a = Wp _1 ' p (£), that is for each given / G L p (0,T; X ) 
there exists a unique solution h G H^(Q, T; X ) n L p (0, T; X\) of the problem 

d t h(t) + A(0)h(t) =/(*), te(0,T), fc(0) = 0, 
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where X x = If U C <" 4/p (S) is a sufficiently small neighborhood 

of zero, then, by a perturbation argument, also A(ho) has maximal L p -regularity, 
whenever ho € U. 

Note that the principal part in (|4.9[) is local in time. Furthermore, by Corollary 
13.41 we have 



Fi, T G C^A^r x B H i nL2 JO,R);L q (0,T;X )), 

for some q > p. This means that the nonlocal term Fi t T is somehow of lower order 
with respect to t. Based on this fact we are in a situation to apply existence and 
uniqueness results for quasilinear evolution equations with main part being local in 
time. We show that the nonlocal term F\ satisfies the Lipschitz estimate 

\\Fi, T {hi) ~ F hT (h 2 )\\ Lp( o,T;X ) < it(T)\\hi - h 2 \\ El T (4.10) 

for all hi, hi £ B r ^ El T , where k(T) — >• 0+ as T — > 0+ and 

B r<ElT := {h e E 1>T : \\h - K\\ ElT < r, h(0) = h }, r £ (0, 1]. 

Here h* £ ^i,t solves the linear Cauchy problem 

d t h*(t)+A(0)h*(t)=0, te{0,T o ), h*(Q) = h , 

for each 7b > 0. Let T £ (0,7b), S > such that ||/io|lyp 4 - 4 /p(£) < ^ follows that 

IIM*)llw^- 4/,> (E) 

< ||ft(t) - /i»(t)|| w 4-4/ P(s) + ||/i*(t) - hoWwf-V'p) + W h °Wwt i/p (v) 

< Mr + sup \\K(t) - h \\ 4-i/ P( . + 5 

teto.ro] Wp w 

< a — e, 

for all /i £ S r EiTi provided that r, 7b, (5 > are sufficiently small. Here a > 
denotes the number in the definition of the set A et R in Theorem 13.21 

Choosing R> r + ||/i*||ei t we obtain that B rtEl T G ^4 e ,_R for all T 6 (0, To). It 
holds that 



Fi, T (/n)-F 1>T (/i 2 ) = 

Hence 



f DF x , T (h 2 + 6{hi 
Jo 



h 2 ))de 



(hi~h 2 ). 



\\F ltT (hx) - F liT (h 2 )\\ Lq{0yT . Xo) < \\Fi t T (e(hi)) - F liTo (e{h 2 ))\\ Lq{0yTo . Xo) 
< CWeihx) - e(h 2 )\\ El , To 



< CM\\hi - h 2 \\ ElT + CM\\h x {Q) - h 2 {Q)\\ 



W p 4 - 4/p (E) 



for all hx,h 2 £ B r ^ 1T , where 

C := SUp{||L> J F 1 i iTo (/l)||£( El . ro ;L 9 (0,To;Xo)) : h G £r,E liT() } > 0. 

and e denotes an appropriate linear extension operator from Ei,t to Ei,t oj T < Tq, 
such that 

\\e(h)\\ EltTo < M (\\h\\ El T + ||/i(0)|| w 4-4 /P(s) ) 

holds for all h £ Ei ; t and M > does not depend on T < 7b and h (see e.g. [6j 
Lemma 7.2]). 
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Since q > p, an application of Holder's inequality yields 

\\Fi, T (hi) - F lyT (h 2 )\\L p (o,T:X (> ) < T pi \\Fi, T (hi) - Fi,t(^2)IU 5 (o,t ; x ) 

<T s ^CM\\h 1 -h 2 \W iT 

for all h\,h2 £ -Br.Ei T - Therefore we can choose n(T) = T 3 ^ CM. In particular, 
the nonlocal term Fi t x(h) is a small perturbation in L p (0,T; Xq) provided that 
T > is small enough. This can be seen as follows 

||*i,tW||mo,T;X ) - W F hT(h) - F hT {K)\\ LA0<T . Xo) + \\F hT (K)\\ Lp{0tT . Xo) 
<K(T)r + T a #\\F 1 , T {h m )\\ Lt[Q , TiXo) , 

for all h £ />', . , T and the right side of the last inequality can be made as small as 
we wish, by decreasing T > 0. 

We may now follow e.g. the lines of the proof of Theorem 2.1] to conclude 
that for each initial value Hq £ U there exists a possibly small T > such that (|4.9[) 
admits a unique solution h £ H^(0,T;X ) nL p (0,T;Xi) which depends (locally) 
Lipschitz continuously on the initial data /iq. 

We have proven the following result. 

Theorem 4.1. Let 3 < p < 2(n + 2)/n, n = 2, 3, R > and U = S T ..4-4/ P ,„, (0, 5). 

W p (2^) 

Then there exist a sufficiently small 5 > and T > smc/j i/iai £/ie (transformed) 
system fl3.6[) - (|3.11[) . (|4. 1 [) - (|4.5[) /ias a unique solution 

(u, q, n, h)EY^xY£x L p (0, T; VK p 2 (0\E)) x Ei, T , 
provided that ho £ ?7 and w £ ^o(^)" n ^2,<r(^), Ikolli/ 1 < -R- 

5. Qualitative Behavior 

This section is devoted to the long-time behavior of solutions to (jl.l[) - (|1.10p 
starting close to equilibria. We will study the spectrum of the full linearization 
of the transformed two-phase Navier-Stokes/Mullins-Sekerka equations around an 
equilibrium. Since, among other things, the divergence-free-condition for the veloc- 
ity field v is destroyed under the Hanzawa transform, we have to split the solutions 
into two parts, one part which is divergence free and the remaining part which is 
not. The treatment of the first part is done by considering the so-called normal 
form of the equations in exponentially weighted spaces and the fact that the set of 
equilibria can be parameterized over the kernel of the linearization. The remaining 
part, which is not divergence free can be handled by the implicit function theorem. 

For simplicity we assume that the dispersive phase is connected. Moreover, we 
assume for simplicity that m = 1. (By a simple scaling in time one can always 
reduce to that case.) Note that the pressure p as well as the chemical potential /i 
may be reconstructed by the semifiow (v(t),T(t)) as follows: 

{Vp\V(j)) L2 = (/i ± Au - v ■ Vv|V<?(>) L9 for all £ W£(n), 

M = n^(Dv)iy m ■ v m j + aH on T(t), 
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and 

mA/i = 0, t>0, i£ ^(t), 

H\ m =<rH, t> 0, x G r(t), 

z/n • mV/j = 0, i > 0, a; G 90. 

Therefore we may concentrate on the set of equilibria £ for the flux (v(i),T(t)) 
which is given by 

£ = {(0, Sr(x q )), Sr(x ) C Q is a sphere}. 

The linearization of the (transformed) two-phase Navier-Stokes-Mullins-Sckcrka 
problem around an equilibrium (0, S) G £ reads as follows: 



d t u - 


p^Au + Vq = /„, 


t > 0, x e fi ± , 






div u = / d , 


t > 0, x e fi ± , 






- a{A^h)vY, = g u , 


t > 0, ieS, 






M = o, 


t > 0, 






u = 0, 


t > 0, xean, 




dth — u ■ 




t > 0, x g s, 


(5.1 




A77 = A, 


t > 0, x g Q*, 






??|s + Ash = g n , 


i > 0, x G S, 






d v r) = 0, 


t > 0, xeffl, 






u(0) = Mo, 


x € Q ± , 






h(0) - ho, 


x e s, 





where .As = 1J j^rI + As and As denotes the Laplace-Beltrami operator on S. 
We want to reformulate (|5.1|l as an abstract evolution equation. To this end we 
introduce the Banach spaces X := L 2 , a (n) x Wp 1/P (E) and X x := {L 2>u (tt) n 
M^ 2 2 (0 \ E) n ) x Wp~ 1/p {i:), where 

L 2)(7 (ft) := {?i G C5°(n) n : div u = 0}"'" i2(ri) . 
Define a linear operator A : -D(A) C X\ — > Xq by means of 

A(u, h) := (-/x ± Au + Vg, -it • - [ft*Tj]), 

with domain 

D(A) = {(u, h) G Xx : u = on 90, [u] = on S}. 

Here g G iJjL (f2\£) and 77 G W^(f2\£) are determined as the solutions of the 
elliptic transmission problems 

(Vg|V0) i2 = (/j^AulVcf)^ for all G H x (n), 

Jg] = 2[ A i ± (i?M)j/ s • j/ s ] + a^ s /i on E, 
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and 

A?7 = 0, i>0,i£l]\S, 

VU + Azh = 0, OO.ieS, 

d v <q = 0, OO.ie dQ. 

In the sequel we will use the solution formula 

Vq = Ti(n ± Au)+T2(2{v ± (Du)vx • M)- 

Setting z = (it, h) and / = (f u ,9h) we may rewrite (15.11) as 

z(t) + Az(t) = f(t), t > 0, z(0) = z := (uo, ho), (5.2) 

provided that fd = g u = /?; = ,9r; = 0. The operator A has the following properties. 

Proposition 5.1. Let n = 2,3, p G (3, 2(n + 2)/n), > 0,<r > 6e constants 
and let Xq and A be defined as above. Then the following assertions hold. 

(1) TTie linear operator —A generates an analytic Co-semigroup e~~ At in Xo 
which has the property of maximal L p -regularity. 

(2) The spectrum of A consists of countably many eigenvalues with finite alge- 
braic multiplicity and is independent of p. 

(3) — A has no eigenvalues A with nonnegative real part other than A = 0. 

(4) X — is a semi-simple eigenvalue with multiplicity n + 1, i.e. Xq = N(A)(B 
R{A). 

(5) The kernel N(A) is isomorphic to the tangent space T Zm £ of £ at the given 
equilibrium z* = (0, E) G £. 

(6) The restriction of e~ At to R(A) is exponentially stable. 

Proof. Consider (JSHJ with f d = g u = /„ = g v = and let J = (0,T), T > 0. 
Suppose that 

h e w${J\ n l p (j- wf-V'pi)) 

is known. Then solve problem (|5.1[l 1 - (j5. ip 5 with initial value uq G (f2) n Hl^cr 
by Theorem I A. 1 1 with g = a = to obtain a unique solution 

u = S hT (h) G ^(0, T; L 2 ^(fi)) n £oo(0, T; n L a (0, T; tf 2 (r! \ E) n ), 

for each T > 0. Plugging it = Si^ih) into (|5.1l) e and denoting by n = S2(—Aj:h) — 
—S^Ash) the unique solution to (|5.1|) 7 8 9 , we obtain the linear nonlocal problem 

d t h - m{d^S 2 {A^h)\ = S hT (h)+g h on E T , h{0) = ho on E. (5.3) 

By [271 Proof of Theorem 4.1] the operator [h >-> m\d vl ,S2{A^h)\] has maximal 
Lp-regularity. Furthermore it holds that Si^Qi) G L q (0,T;Wl ^ P (E)) for some 
q > p which means that this term is of lower order in L p (0,T; Xq) compared to 
mld Us S2(Ash)J\. This can be seen as in the proof of Corollary 13.41 

Hence, by perturbation arguments we may conclude that (|5.3p has for each given 
ho G Wp~ 4 ^ p (E) a unique solution 

h G W}(J; Wp _1/,p (E)) n L P (J; ^"^(E)). 

In other words, we have shown that for each T > and for each given / = (/„, gh) G 
L%{J; L2,<r(^)) x L P (J; Wp ^ P (E)) there exists a unique solution z — (u, h) of (|5.2j) 
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with 

U G x (0, T; L 2 , CT (ft)) R Loo(0, T; H^(fl) n ) n L 2 (Q, T; ff 2 (fl \ E)™), 

and 

provided that u a G £ 2 ,<r(^) n F^(fi)" and h £ Wp~ 4/p (S). 

Mimicking the proof of [551 Proposition 1.2] it follows that the operator —A 
generates an analytic semigroup in Xq = L 2)CT (f2) x Wp 1//p (E). 

By compact embedding, the resolvent of A is compact and therefore the spectrum 
<j{A) of A consists of countably many eigenvalues with finite algebraic multiplic- 
ity and a(A) does not depend on p, by classical results. Let A G er(— A) with 
eigenfunctions (u, h). Then the corresponding eigenvalue problem is 

Ait — Att + Vg 
div u 

-21^(^)1^ + [g]i*3 - <r(„4 s % s 

M 

A/i-u-^-I^??] =0, zeE, (r,.4) 
J?|s + -4u/» 

Taking the inner product of equation ()5.4p i with it, integrating by parts and invok- 
ing the boundary as well as the transmission conditions, wc obtain 

A||«||l + n^DuWl + a\\Vr)g - a\(Avh\h) L2(s) = 0. (5.5) 

If A ^ 0, then 

A / hdo — / vs ■ udo + / \d v ^rj\do = I divudx = 0, 
Js Js Js Jn+ 

hence h has mean value zero. It is well-known that the operator As = ^ + As is 
negative definite on L 2 ,(o)(^)- Taking real parts in (|5.5j) it follows that ?y = const 
and -Du = 0, hence u = by Korn's inequality since it|an = 0. This in turn yields 
h = by (|5.4p £ . showing that there are no eigenvalues A ^ of —A with Re A > 0. 
Next we show that A = is an eigenvalue of A. If A = 0, then (|5.5I) implies 
n = rjoo = const and Du — 0. Hence, as before, u — by Korn's inequality. Since 
q is constant by (|5.4p i it follows from (|5.4[U 7 that 

r]oo = — = —5n~h + Ash, 

which is a linear second order partial differential equation for on E. Note that 
a special solution to this linear equation is given by the constant function = 
r/ooR /(n — 1). The solution space C of the corresponding homogeneous equation 
Ash = is given by 

C = span{Yi, . . . ,Y„}, 



0. 


X 


e n\E, 


0. 


X 


g n\s, 


0. 


X 


gE, 


0. 


X 


GE, 


0. 


X 


GS, 


0. 


X 


g q\e, 


0. 


X 


GE, 


0. 


X 


g on, 


0. 


X 


e an. 
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where Yj, j G {1, . . . , n}, are the spherical harmonics of degree one. Furthermore 
it holds that dim£ = n. Since the constant 7700 = M/o - is arbitrary, we see that 
dim N(A) = n + 1. 

Let z\ G N(A) such that Az = z\. The corresponding problem for z — (u, h) T is 
given by 





~Au + Vq 


= 0, 






div it 


-0, 


x g \ £, 


2\jiDu]vE + 


lq]i/ s - a{AY,h)v T , 


-0, 


ieS, 




H 


= 0, 


ieS, 






--h u 


xeS, 




A77 


= 0, 


x e Q\E, 




»7 |s + Ash 


-0, 






d u T] 


= 0, 


a; G 90, 




u 


= 0, 


x g an, 



(5.6) 



since Z\ = (0, /ii) and hi — X)J=o ^0 := F rom t ne divergence condition we 

obtain 

= / div -u eta; = / (m • i/£ + [<9„ E 77] ) do = — / hido, 
■I" -h: ■>'■: 

hence hi has mean value zero and this in turn implies A-^hi = 0. Multiplying (|5.6p i 

by u, integrating by parts and taking into account the boundary and transmission 

conditions, we obtain 

2\\^Du\\ 2 2 + a\\Vr)\\l + <t(-4 s ^i)l 2 (£) - 0. (5.7) 

Since Ay, is self-adjoint in L2CE) it follows that the last term in (|5.7|) vanishes and 
then, as before, 77 = const and u = 0, by Korn's inequality. In this case (j5.6|) ; 
yields hi = 0, i.e. z G N(A), hence N(A 2 ) = N(A). Since A has compact resolvent, 
it follows that R(A) is closed in Xq and A = is a pole of (A — A)^ 1 . Therefore 
p51 Remark A. 2. 4] yields that A = is semisimple, in particular it holds that 
Xq = N(A) © R(A), Moreover, the restricted semigroup e~ At \n(A) is exponentially 
stable, since we have a spectral gap. 

Finally we show that the tangent space T Zt £ of £ at z* = (0, E) G £ coincides 
with N(A). This can be seen as follows. Assume w.l.o.g. that E is centered at the 
origin of E™ with radius R. Suppose S is a sphere that is sufficiently close to E. 
Denote by (yi , . . . , y n ) the center of S and let R + yo be the corresponding radius of 
S. Then by |12[ Section 6] the sphere S can be parametrized over E by the distance 
function 



\ 



5>r, +(i?+yo) 2 -£yf. 

vi=i / 3=1 



Denoting by O a sufficiently small neighborhood of in the mapping d : 

O — > Wp~ l ^ p (Ti) is smooth and the derivative at is given by 

n 

d'(0)w = WjYj, for all w G M™ +1 . (5.8) 
3=0 
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Therefore, near £, the set of equilibria £ is a smooth manifold in X\ of dimension 
n + 1 and T Zr £ = N(A) by (|5T5j) . 

Since X = N(A) © R(A) and cr(A| fl ( A) ) C C+ it follows that the restricted 
semigroup e~ \ma) 1S exponentially stable. The proof is complete. 

□ 

We are now ready to prove the main result of this section. Note that the trans- 
formed equations near an equilibrium (0, X) 6 £ read as follows. 



8 t U — [J, Au + V-7T = F u (u, 7T, h), 


t > o, x e n ± , 


div u — Fd(u, /i), 


t>o, x e fi ± , 


[/x(Vu + Vm t )]^ + [tt]i/£ - a(Ash)vT, = G u (u, h), 


i > 0, x e S, 


M - o, 


t > 0, x G S, 


dth-u-VT, - [<9„ E ?/] = Gh(u,r), h), 


i > 0, 


Ar) = F v (r],h), 


t > 0, is!! 1 , (5.9) 


ri\v + Axh = Gr,(h), 


t > 0, x g S, 


d v n = 0, 


* > o, x e an, 


u = 0, 


* > o, x g an, 


tt(0) = M , 


x g n ± , 


MO) = fto, 


x G S, 



where the derivatives of the nonlinearities on the right hand side with respect to 
(u, h) vanish at (u, h) = (0, 0) for constant tt and constant r\. 

Theorem 5.2. The equilibrium (0, X) G £ is stable in the sense that for each e > 
there exists some 5(e) > such that for all initial values (uo,ho) subject to 

||/i || w *-4/ P(E) + IKHffi(Q)" < 5(e) 

there exists a unique global solution (u(t),h(t)) of (|5.9p and it satisfies 

WHtfWw*-*/*^ + \\u(t)\\ H i isl)n < e for all t > 0. 

Moreover, there exists some hoc G Wp such that O/j^E = dBn(x) C O, for 

some i? > 0, x G n, and 

/i™ (iW*) - ft °°Hw*- 4/f, (S) + IK*)llHj(n)») = °- 
T7ie convergence is at an exponential rate. 

Proof. The nonlinear phase manifold for the semiflow is given by 

VM = {(u,h) G H%(p,) n x Wp 4 - 4/p (E) : divu = F d (u,h)}. 
In a first step we want to parametrize VAA over its tangent space at (0,0), that is 

VM := {(u, h) G #o(fi)" x % 4 ~ 4/p (X) : divu = 0}. 
To this end we consider the generalized Stokes equation 

-Am + Vtt = in n, 

div« = / inn, (5.10) 

u = o on an 
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for which we have the following existence and uniqueness result. 
Proposition 5.3. For every f G iJLJVt) = {u G L 2 (Q) : J n udx = 0} the Stokes 
problem (|5.10[) admits a unique solution (u,ir) G x -^(o)(^)> which depends 

continuously on f G L^(f2). 

Proof. The proposition is a special case of [331 Theorem III. 1.4.1]. □ 

With the help of this result we may continue as follows. For a given (u, ft) G VMo 
with a sufficiently small norm, we solve the auxiliary problem 

-Au + Vtt = in fl, 

divu = P F d (u + u,h) infi, (5.11) 

u = on 9fi, 

where Pq : L 2 (fl) — > i^ ^ (SI) is defined by PqJ = f — j^j Jq fdx. Since the Frechet 
derivatives of the nonlinearities vanish in (0, 0), the implicit function theorem yields 
the existence of a ball -B(0, r) C Hq (fl) n D Wp 4 ^ P (E) and a unique solution 

(u,7f) = G i?o( fi ) x L (0)(°) 

with a function G C^fl^r)) such that 0'(O) = 0. Define 9^ (a;) as in (l2~8|) with 
ft. replaced by ft, which does not depend on t. Let v(x) := (u + u)(Q^ (x)). Then 
v G i^(O)" and 

divv(x) = Tr[D6^ T (x)V(u + w)(e^ 1 (a;))] 

= Tr[(L>6^ T (x) - I)W(u + u)(9r 1 ( i r))] + div(w + u)(9^ 1 (x)) 

= -F d (u + u, h)(O h \x)) + P F d (u + u, h)(G^(x)) 



since divu = 0. Because of = J n divv(x)dx, it follows that PoFd(u +u,h) = 
F d (u + u, h). 
Let 

P: ff^O)" x Lf 0) (O) -> ^(O)", P(u,tt) = u, 

and set (/>(«,, ft) = P<j>(u, ft). It is not difficult to see that 

(u, ft) := (u, ft) + (</»(u, ft), 0) G VM. 

Note that this mapping is injective. For the final construction of the parametriza- 
tion we have to show that this mapping is also surjective. For that purpose we 
solve the linear problem 

-Au + Vtt = mfi, 
divu = P F d (u,h) in Q, (5.12) 
u — on dfl, 

for given functions (u, ft) G VM.. Setting [u, ft) = (u — u, ft) we obtain that u G 
ff^(fi) n and 

divu = F d (u, ft) - P F d (u, ft) = y F<j(u, ft)dx. 
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Since = J n div udx this yields PoFd(u, h) — Fd{u, h). 

Furthermore it holds that (u, h) G PAio and u = 4>{u, h) by injectivity. This 
in turn proves surjectivity. Observe that also <j>(0) = 0. This can be seen as 
follows. Suppose that u = h = 0. Then obviously u = and ff = const, is a 
solution of (|5.1ip . By the uniqueness it follows that (f)(0) = 0. Furthermore, if 
) is an equilibrium of (15.91) , then = and 

Voo = [tTooI/ct = H(hoo) = const. 

Since F^O,/^) = 0, the unique solvability of (|5.1ip implies that (f>(0, h x ) = 0. 
This is reasonable since the equilibria are contained in the linear phase manifold 
VM - 

Let (uo,ho) = (uo,ho) + (<p(uo,ho),0) G VM. and let (u,h,Tr,rj) be the solution 
of (|5.9p to this initial value on some interval [0, a]. With the help of the map <fi we 
want to derive a decomposition for (u, h). To be precise we want to write 

(u, h) = (u x , hoo) + (u, h) + (u, h), 

where (u,h)(t) G VMq for all t G [0, a] and (u^, hoo, Tco, Voo) is an equilibrium of 
(|5.9[) . Consider the two coupled systems 

uiu + dtu — n Au + V7f = F u (u, 7T, h) 

div u = Fd(u, h) 

-P S MV?1 + Vu T )K = G T (u,&) 

-([^(Vu + Vu T )K|^) + [tt] - oA^h = G„{u, h) + G 7 (h) - G 7 (/ioo) 

[«]=Q 

u\an = (5.13) 
ujh + d t h -u-vs- [d Vs ffl = Gh(u,n, h) 
Afj = Fr,{r],h) 
fj\s + Avh = G v {h) - G v {hoo) 
duf)\dn = 
6(0) =(/>(uo,h ), h{0)=0, 

and 

d t u — /i ± Au + V7r 

div -it 

-P s \H ± {Vu + Vu t )K 

-(MVu + Vu T )p s ) ■ ve + [tt] - crA s h 

N 
u 

d t h-u-vx- Id^fj} 
Afj 

fj\ s + Avh 



LOU 


t > 


o, 


x e fi\E 





t > 


o, 


x e n\s 





t > 


o, 


x e E 





t > 


o, 







t > 


o, 


ieE 





t > 


o, 




u)h 


t > 


o, 


ieS 





t > 


o, 







t > 


o, 


ieE 





t > 


0. 


x g 9£7 



(5.14) 
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with initial values u(0) — uo and h(0) = ho — hoo. Here 77 = ix x + 77 + 77 and 
i] = ?7oo + fj + t)- We recall that Uoo = and 7700 , 7700 are constants and it holds that 

??oo = [tTooJ/ct = H(hoo). 

With the help of the operator A introduced above, we may rewrite problem 
(|5+4| as 

z(t) + Az(t) = R(z)(t) for t g (0, T), z(0) = z - z^, (5.15) 

where zq :— (uo,ho), z x :— (0,/ioo), z = (u,h) and z = (u,h). Here the mapping 
R is given by 

R(z) = (u(I -Ti)u,uh). 

Thanks to Proposition 15 . 1 1 we have the decomposition Xq = N(A) ® R(A). Let P c 
denote the spectral projection corresponding to a c (A) — {0} and set P s = I — P c . 
Then R(P C ) = N(A) and R(P S ) = R(A). Following the lines of [23 Remark 
2.2 (b)] we may parametrize the set of equilibria near over N(A) via a C 2 map 
[x i-> x + ip(x)] such that ^(0) = tp'(0) = and R(ip) C D(A a ), where A s = AP S . 
This is true, since the nonlinearities on the right side in (|5.9I) are bilinear and 
smooth. 

For sufficiently close to there exists x^ such that z M := x^ + ^(x^). 
Introducing the new variables x := P c z and 

y := P S S - iPixoo + P c z) + V>(xoo) 

we obtain from (]5 . 15[) the so-called normal form 

k = T(z), x(0)=x -x oo , 

y + A s y = S(xoo,x,z), y(0) =y , (5.16) 

where T(z) = P c R(z), 

5( Xoo ,x, z) = P s R{z) - ^(^(xoo + x) - ^( Xoo )) - ^'(xoo + x)T(z), 

and x := P C S , yo := P s S n — tpfa) with z — (uo,ho). Observe that 5(0) = 
5'(0) = by the properties of the function ijj and since T(0) = 0. 
Let 

Ei(M+) := H 1 (R+;L 2 .a(n)) nL 2 (I + ;ff 2 (fl\E)), 
E 2 (K+) := W^(K+; W^(E)) n L P (R+; <" 1/p (S)). 

and let 

Ei(K+,5) :={neL 2 (R + ;i2(fi)) : e 5t w GEi(R+)}, 

E 2 (K+,tf) := {v g L p (R+;L p (fi)) : A g E 2 (M+)}, 

where <5 g (0, Jo) and <5o > depends on the spectral bound on the operator A s 
(see Proposition El]). Clearly, T : E(R+,5) -> E c (R+,<5), where 

E(R+ , 5) := Ei (R+ , S) X E 2 (R+ , 5) 

and E c (R+,<5) := P c E(R+,(5). For given (x ,y ,z) we want to solve (|5~T6| for 
(x, y, Xqo). First, for given (x ,z) £ X^ x E(R + ,<5) with := F C X we define 

/>oo 

xoo := x + / T(z)(s)ds =: Kx(x , z) g X C . 
Jo 

Then 

x(i) := - / T{z){s)ds =: K 2 {z) 
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solves the first differential equation in (|5.16[) and 

r-OO 

x(0) = - / T(z)(s)ds = x Q -x 00 . 
Jo 

Observe that by Young's inequality we have 

x e P c [H 1 {R+,6;L 2 An)) x W p 1 (E+, 5; W^ P (T))]. 

These exponentially weighted function spaces are defined in exactly the same way 
as Ej(R+,5). 

Substituting the expressions for Xoo and x into the function 5, we obtain 
y + A s y = Si(x ,z), y(0) = y , 

where 

Si(xq,z) := S(K! (x ,z) 7 K 2 (z),z) , 

and yo € Xq n VA4q. If one takes into account that the first component of -0 is 
identically zero, it follows easily from the definition of S and the smoothness of ip 
that 

Si(x 0l z) : X C Q x E(K+,$) -> X s (M + ,<5) 

where 

X fl (K + ,«5) :=P s [L 2 (R + ,6;L 2 , a (n)) x L p (R + ,S;W^ 1 / p (^))]. 
Here Xq := P S X$. Since a(A s ) C C+, we obtain for 5 > sufficiently small 

y= f^ + A s ,tA (5i(xo,«),y ) eE s (R+,<5), 
where tr v := v(0) and 

E s (M+,(5) := P s E(R+,5). 

Here <5 > depends on the growth bound of the semigroup. Putting things together, 
we see that 

z = <?(x , y , z) := x + ip(x + x^) - ^(xoo) + y 

and 

Zoo = Goo(Xo,yo,z) := Xoo +^(X 00 ). 

We turn our attention to (|5 . 13|) . Let L w be the linear operator defined by the left 
side of (|5.13|) . Then we can rewrite (|5.13j) in the shorter form 

L u w = N(w oa + w + w) - N(w oa ), (5.17) 

with initial value w(0) = wq '■= (4>(uo,ho),0). 

Here we have set Woo — (uoo, hoc, 0, 0). Due to the first part of the proof, the 
nonlinearities on the right hand side of (|5.13l) depend only on (x ,yo,w), where 
to = (u, h, 7f, fj) since = (C?oo( x o, Yo, h), 0, 0) and since there exists a function 
H such that 

to = (u, h, ft, fj) = H(x , y , u, h). 
This follows from the considerations above, as (ft, fj) can be written in terms of 
(u, h) and (u, h) = z = G(xq, yo, z). Moreover, the right hand sides in (|5.13|) do not 
depend on (iToo, V<x), since these quantities are constant. 
In order to solve (|5.13l) we define 

M (x , y , to) := N( Woo +w + ext s [(<p(u Q , ho), 0) - (6(0), h(0))] + to) - N(w oa ), 
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where 

ext 5 : {^(SJfnivP)} x % 4 " 4 /f(E) -> E(R+,S) x {0} x {0}, 
such that (extsz)(0) = (zi, Zi, 0, 0) with z = (zi, z%). We want to solve the equation 
L w w = M(x ,y ,u>), (wi,w 2 )(0) = ((f)(ua,ha),0), (5.18) 
by the implicit function theorem. Let 

E(R+,6) :=E(R+,<J) x L 2 (R+, 5; ff^^E)) x L p (R+, 5; VF p 2 (0\S)) 
and define 

if(x ,y ,w) := u) - (L w , tr) _1 (Af(x , y , w), (<t>(u , ho), 0)). 

The mapping X : B(r, (5) — > E(R+,<5) is well defined, provided that u > is suf- 
ficiently large since (M(xo, yo, (</>(tto, ^o); 0)) satisfies all relevant compatibility 
conditions at t = 0. To be precise, we have [</>(ito, ho)} — 0, 0(uo, /io)|an = as 
well as 

div0(u o , ho) = div(u + 4>{u , ho)) = F d (u + 4>(u , ho), ho), 
since div uq = 0. Here we have set 

M(r,S) := 

{(xo,y ,«)) € A C x Wn?M ) x E(M+,<J) : ||(x ,y , ^II^^e^) < r}, 

where r > is sufficiently small. 

Note that M (0, 0, 0) = M'(0, 0, 0) = since 0(0, 0) = 0'(O, 0) = 0. Therefore the 
implicit function theorem yields a ball 

B(0,rfcl c x (X s nPA4 ) 

and a unique solution w = <£>(x ,yo) of (|5.18[> . where $ 6 C 1 (B(0, p)). Note that 
by construction, w is a solution of (|5.13l) . 

Finally this shows that (u(t),h(t)) as well as (u(t),h(t)) converge in _ffg(il) n x 
<" 4/p (E) to zero as t tends to infinity at an exponential rate. 

Therefore (u(t), h(t)) -> (wqo, /loo) in H^(Q) n x Wp~ 4 ' p (S) as t -> oo, where the 
equilibrium (itoo, /loo) is determined by (uq, ho). □ 

Appendix A. Maximal Regularity for the linear Stokes System 

For the following let Qg , be bounded domains with C 3 -boundary such that 
r := dfl+ C n and let fl = Q \ ft+ . Recall that H^ Q) (Q) = H 1 ^) n i 2 ,(o)(^) 
and J ff (0 J(fi)=^ 0) (O). 

In this appendix we consider the unique solvability of the system 



d t u- /i ± Au + Vg = / in ^ x (0,T), (A.l) 

divu = .g in fi^ x (0,T), (A.2) 

M=Q onr x(0,T), (A.3) 

[vr a ■ T(u, q)j = a on T x (0, T) =: T ,t, (A.4) 

u|an o =0 onan x(0,T), (A.5) 

u|i=o = «o on n , (A. 6) 
where T(u, q) — ^Du — ql in Qq. 
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Theorem A.l. Let < T < Tq < oo, n > 2, and fl C R ra be a bounded domain 
with C z -boundary. For every v <= H&(tt) n ,f € L 2 (Q T ) n 1 g £ L 2 (0, T; H 1 ^ \ Tq)) 
withg&H l {0,T-H-^(n)), 

a G Hi(0,T;L 2 (T ))nL 2 (0,T;Hi{T )) =: H^(T , T ) 

such that 

div t>o = g\t=o, / g{t,x)dx = for almost all t G (0, T) (A. 7) 

i/ie system (|A.1[) - (|A.6[) /ias a unique solution 

u G ff^O, T; L 2 (ft) n ) n Loo(0, T; ^(fi)") D L 2 (0, T; 2 (ft \ r )«) 

Moreover, there is some constant C independent of T G (0, To], w, f,g,a,vo such 
that 

||(a tU) v u )|u 3(J;i3( n ) ) + X)ll( v Vv?)|| i2(J . i2(n ± )) 

± 

< C g (||(/,Vg)|| i2( g r) + ||5 t <?|| L2(J . ff -i (n)) + ||a|| H i 4(ro ^ + |M|jp ( n)) 
where J = [0, T] . 

Remark A. 2. The result follows from a result announced by Shimizu [31], where 
a general L 9 -theory is discussed. In the case q — 2, the proof is much simpler 
since Hilbert-space methods are available and the result basically follows from the 
resolvent estimate proved by Shibata and Shimizu in [30] . For the convenience of 
the reader we include a proof. 

Proof of Theorem IA.lt First we consider the case g = a = vq = 0. We can 

assume without loss of generality that / G £2(0, T\ L 2 ,<r(^))- Otherwise we replace 
/ by P a f and q by q — qi, where Vqi = f — P a f . Then (|A.1[) - (IA.6I) are equivalent 
to the abstract evolution equation 

ju{t)+Au{t) = f(t), te(0,T), (A.8) 
u(0) = 0, (A.9) 

where A: V{A) -> Lz t<T {&) with 
Am|q± = — v^Au + 'Vq 
27(A) = (u G ffo^ )" n L 2, CT (0) : V 2 u| ± G L 2 (n±), [2v ■ ^Dv] T = o} 

where q G L 2 ,(o)(^) with Vq| f2 ± G £2(^0)™ is uniquely defined by 

Aq = in 
M = l^d^j on r , 

d v q\on = v ■ fi~Au\dn on <9fi. 

Because of [30l Theorem 1.1], A is a generator of an exponentially decaying analytic 
semi-group and the graph norm ||it||x>(A) is equivalent to 

\\u\\ H i(n) + J2\\ v2u h 2 (nty 
± 
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Since L2,<r(^) is a Hilbert space, for every / 6 £2(0, T; £2,0-) there is a unique 
u e H 1 {0,T;L 2 ,a) nL 2 (0,T;D(A)) solving (|X8|) -([A ^ and 



eft 



+ ||Au 

i2(0,T;L 2 ,„) 



li2(0,T;i 2 , CT ) ^ C||/l|i2(0,T;i, 2 , CT ) 



with some C > independent of T £ (0, 00]. In the case T = 00 this statement 
follows from [9] or (TTJ Theorem 4.4], part "(ii) implies (i)", where we note that 
7?.-boundedness of an operator family on a Hilbert space coincides with uniform 
boundedness, cf. [TTJ Section 3.1]. The case T < 00 follows from the latter case by 
extending / : (0, T) — > H by zero to / : (0, 00) — > H . This proves the theorem in 
the case g = a = vq = 0. 

The general case can be reduced to the latter case as follows: First we reduce to 
the case (/, 5,«o)ln+ =0. To this end let 

v+ e H\0,T;L 2 (n+) n )nL 2 (0,T;H 2 (n+r), q+ e L 2 (0, T; H 

be the solution of 

d t v+ - fi + Av + + V<7+ = /| n + in n+ x (0, T), 

divv+=g\ nS infi+x(0,T), 

v To • (2fi+Dv+ -q+)=0 on T x (0, T), 

v + \ t=0 =v \ n + infi+. 

The existence of such a v + follows from well known results for the instationary 
Stokes system with Neumann boundary conditions, cf. e.g. [3]. Moreover, there is 
some constant C > such that for every < T < 00 

\\(d t v+,Vv+,^v+,Vq+)\\ L2{J . L2{ni)) 

< Cg (lllf, Vff)|| ia(Jxn +) + ||fltfllL a (j ;ff -i(n+)) + \\ v o\\m(n+)) ■ 
Now we extend v + and q + to some functions 

v + €La(0,T;H a (Q o )) n nH 1 {0,T-,La(Sl o )) n , q+ e i 2 (0, T; H\Q+)) 

satisfying an analoguous estimate as before. Now subtracting (v + ,q + ) from (u,q) 
we reduce to the case (f,g,vo)\ n + = 0. Next we observe that 

g\ n - GH'^T-H^Q-)) 

because of 

g(x,t)ip(x) dx = I g(x,i)(p(x) dx 





for every ip S HhJD,^ ), where (p s iJjL(fio) is an arbitrary extension of <£> to f2o- 
Hence there are some 

w-GH 1 (0,T;L 2 (fio) B )nL 2 (0,T;H 2 (no) n ) ) <T E £ 2 (0, T; H\% )) 
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solving 

d t v- - fx-Av- + Vq~ = /| n - in % x (0,T), 

divv- = g\ n - inft o x(0,T), 

v~ = on r x(0,T), 

v~\ t =o = v Q \ n - inO . 

Existence of such a solution together with analoguous estimates as for (v + ,q + ) 
follows e.g. from [3j [T3J [16l [TTl [34] . Now extending v~ and q~ by zero to and 
subtracting the extensions from (u,p) we can reduce to the case (/, <?, «o) = 0. 

In order to reduce to the case, where also a T = 0, we construct some A G 
H 1 ^, T; L 2 (n+)) n L 2 (0, T; H 2 {n+)) such that 

|[(A9tAVAV 2 A)|| r < C\\a\\ ii 

IIV ' 1 ' ' J "L 2 (Q T ) — II Hjj4'2 (r ,T) 

and 

A\ t =o = A\r = 0, (^ r „ •2^+DA) T |r = a T) divA = in Q . 
This can be done as follows: Choose some 

A e ff x (0, T; i 2 (0+)") D L a (0, T; ff 2 (f#) n ) 

such that 

||(l,9 t I,Vl,V 2 i)|| i2( n o+x(0 , r)) < C||a|| Hi , i(roT) 

and 

A| ro = i| t=0 = 0, Or • 2M + -Dl) T |r = o T , divl|r = 0. 

The existence of such an A e.g. follows from [3J Lemma 2.4]. Moreover, C > in 
the estimate above can be chosen independently of < T < To for any To > 0. 
Since div A\r = 0, div A G Hq(Qq) and we can apply the Bogovski operator B, cf. 
e.g. [14 , to div I. Then we obtain S(divl) 6 L 2 (J; H 2 (n+) n L 2 {Q) {n+)) and 

H- B ( di V^)|lz, 2 (J;H=(n+)) < C 'll^llL 2 (J;^(f2+))- 

Moreover, due to [15j Theorem 2.5] we also have 

||-B(divI)|| H i (J . M n+ )) < C\\divA\\ H1{J . H -i {n + )} < C'\\A\\ H1{J . L2{n + )y 

Since the Bogovski operator is independent of time, the latter constant can be 
chosen independently of < T < Tq for any To > 0. Altogether, we obtain that 
A := A — £?(div A) has the properties stated above. Replacing nby«- ^Xn + > we 
can finally reduce to the case vq = g = a T = 0. Finally, we can also reduce to the 
case a v = by substracting a suitable extension of a„ from the pressure q. ■ 

Acknowledgments: We are grateful to the anonymous referees and Stefan Schau- 
beck for careful reading previous versions of this work and many comments, which 
improved the paper. Moreover, the authors acknowledge support from the German 
Science Foundation through Grant Nos. AB285/3-1 and AB285/4-1. 

M.W. would like to express his thanks to Gieri Simonett for inspiring discussions 
concerning the proof of the stability result. 



TWO-PHASE NAVIER-STOKES-MULLINS-SEKERKA EQUATIONS 



33 



References 

1. H. Abels, The initial value problem for the Navier-Stokes equations with a free surface in 
Li-Sobolev spaces, Adv. Diff. Eq. Vol. 10, No. 1 (2005), 45-64. 

2. H. Abels, On a diffuse interface model for two-phase flows of viscous, incompressible fluids 
with matched densities, Arch. Rat. Mech. Anal. 194 (2009), no. 2, 463-506. 

3. , Nonstationary Stokes system with variable viscosity in bounded and unbounded do- 
mains, Discrete Contin. Dyn. Syst. Scr. S 3 (2010), no. 2, 141-157. MR MR2610556 

4. H. Abels, H. Garcke, and G. Griin, Thermodynamically consistent, frame indifferent dif- 
fuse interface models for incompressible two-phase flows with different densities, Preprint, 
arXiv:/1104.1336, to appear in "Math. Models Methods Appl. Sci." (2011). 

5. H. Abels and M. Roger, Existence of weak solutions for a non-classical sharp interface model 
for a two-phase flow of viscous, incompressible fluids, Ann. Inst. H. Poincare Anal. Non 
Lineaire 26 (2009), 2403-2424. 

6. H. Amann, Quasilinear parabolic problems via maximal regularity, Adv. Differential Equations 
10 (2005), no. 10, 1081-1110. MR MR2162362 (2006d:34125) 

7. J. Bergh and J. Lofstrom, Interpolation spaces, Springer, Berlin - Heidelberg - New York, 
1976. 

8. F. Boyer, Mathematical study of multi-phase flow under shear through order parameter for- 
mulation, Asymptot. Anal. 20 (1999), no. 2, 175-212. MR MR1700669 (2000g:35166) 

9. L. de Simon, Un'applicazione della teoria degli integrali singolari alio studio delle equazioni 
differenziali lineari astratte del primo ordine, Rend. Sem. Mat. Univ. Padova 34 (1964), 205- 
223. MR 0176192 (31 #467) 

10. I. V. Denisova and V. A. Solonnikov, Solvability in Holder spaces of a model initial- boundary 
value problem generated by a problem on the motion of two fluids, Zap. Nauchn. Sem. 
Leningrad. Otdel. Mat. Inst. Steklov. (LOMI) 188 (1991), no. Kraev. Zadachi Mat. Fiz. i 
Smezh. Voprosy Tcor. Funktsii. 22, 5-44, 186. MR MR1111467 (92c:35126) 

11. R. Denk, M. Hieber, and J. Priiss, TZ-boundedness, Fourier multipliers and problems of elliptic 
and parabolic type, Mem. Amcr. Math. Soc. 166 (2003), no. 788, viii+114. MR MR2006641 
(2004i:35002) 

12. J. Escher and G. Simonctt, A center manifold analysis for the Mullins-Sekerka model, J. 
Differential Equations 143 (1998), no. 2, 267-292. MR MR1607952 (98m:35228) 

13. A. Frohlich, The Stokes operator in weighted L q -spaces. II. Weighted resolvent estimates and 
maximal L v -regularity, Math. Ann. 339 (2007), no. 2, 287-316. MR MR2324721 (2008i:35188) 

14. G. P. Galdi, An introduction to the mathematical theory of the Navier-Stokes equations, 
volume 1, Springer, Berlin - Heidelberg - New York, 1994. 

15. M. GciBert, H. Heck, and M. Hieber, On the equation divu = g and Bogovskii's operator 
in Sobolev spaces of negative order, Partial differential equations and functional analysis, 
Oper. Theory Adv. Appl., vol. 168, Birkhauser, Basel, 2006, pp. 113-121. MR MR2240056 
(2007k: 35034) 

16. M. Geissert, M. Hess, M. Hieber, C. Schwarz, and K. Stavrakidis, Maximal L p — L q -estimates 
for the Stokes equation: a short proof of Solonnikov's theorem, J. Math. Fluid Mech. 12 
(2010), no. 1, 47-60. MR MR2602914 

17. Y. Giga and H. Sohr, Abstract IP estimates for the cauchy problem with applications to the 
navier-stokes equations in exterior domains, J. Funct. Anal. 102 (1991), 72—94. 

18. M. Grasselli and D. Prazak, Longtime behavior of a diffuse interface model for binary fluid 
mixtures with shear dependent viscosity, Interfaces Free Bound. 13 (2011), no. 4, 507-530. 
MR 2863469 

19. M. E. Gurtin, D. Polignone, and J. Vifials, Two-phase binary fluids and immiscible fluids 
described by an order parameter, Math. Models Methods Appl. Sci. 6 (1996), no. 6, 815-831. 
MR MR1404829 (99e:76123) 

20. P.C. Hohcnbcrg and B.I. Halperin, Theory of dynamic critical phenomena., Rev. Mod. Phys. 
49 (1977), 435-479. 

21. J. Johnsen, Pointwise multiplication of Besov and Triebel-Lizorkin spaces, Math. Nachr. 175 
(1995), 85-133. MR MR1355014 (97a:46042) 

22. N. Kim, L. Consiglieri, and J. F. Rodrigucs, On non-Newtonian incompressible fluids with 
phase transitions, Math. Methods Appl. Sci. 29 (2006), no. 13, 1523-1541. MR MR2249576 
(2007g:76008) 



34 



HELMUT ABELS AND MATHIAS WILKE 



23. M. Kohne, J. Priiss, and M. Wilke, Qualitative behaviour of solutions for the two-phase 
Navier-Stokes equations with surface tension, to appear in Math. Ann. 

24. , On quasilinear parabolic evolution equations in weighted L p -spaces, J. Evol. Equ. 10 

(2010), no. 2, 443-463. MR 2643804 

25. A. Lunardi, Analytic semigroups and optimal regularity in parabolic problems, Birkhauscr, 
Basel - Boston - Berlin, 1995. 

26. J. Priiss, Maximal regularity for evolution equations in L p -spaces, Conf. Semin. Mat. Univ. 
Bari (2002), no. 285, 1-39 (2003). MR 1988408 (2004k:35232) 

27. J. Priiss, G. Simonett, and R. Zacher, On convergence of solutions to equilibria for quasilinear 
parabolic problems, J. Differential Equations 246 (2009), no. 10, 3902-3931. MR 2514730 
(2010d:34123) 

28. I. Rousar and E.B. Naumann, Spinodal decomposition with surface tension driven flows, 
Chemical Engineering Communications 105 (2010), no. 1, 77-98. 

29. T. Runst, Mapping properties of nonlinear operators in spaces of Triebel-Lizorkin and Besov 
type, Anal. Math. 12 (1986), no. 4, 313-346. MR MR877164 (88f:46079) 

30. Y. Shibata and S. Shimizu, On a resolvent estimate of the interface problem for the 
Stokes system in a bounded domain, J. Differential Equations 191 (2003), no. 2, 408-444. 
MR MR1978384 (2004d:35203) 

31. S. Shimizu, Maximal regularity and viscous incompressible flows with free interface, Parabolic 
and Navier-Stokes equations. Part 2, Banach Center Publ., vol. 81, Polish Acad. Sci. Inst. 
Math., Warsaw, 2008, pp. 471-480. MR MR2548875 

32. J. Simon, Sobolev, Besov and Nikol'skii fractional spaces: imbeddings and comparisons 
for vector valued spaces on an interval, Ann. Mat. Pura Appl. (4) 157 (1990), 117—148. 
MR 1108473 (92c:46075) 

33. H. Sohr, The Navier-Stokes equations, Birkhauser Advanced Texts: Basler Lehrbiichcr. 
[Birkhauscr Advanced Texts: Basel Textbooks], Birkhauser Verlag, Basel, 2001, An elemen- 
tary functional analytic approach. MR MR1928881 (2004b:35265) 

34. V. A. Solonnikov, Estimates for solutions of a non- stationary linearized system of Navier- 
Stokes equations, Trudy Mat. Inst. Steklov. 70 (1964), 213-317. MR MR0171094 (30 #1325) 

35. E. M. Stein, Singular integrals and differentiability properties of functions, Princeton Hall 
Press, Princeton, New Jersey, 1970. 

36. H. Triebel, Interpolation theory, function spaces, differential operators, North-Holland Pub- 
lishing Company, Amsterdam, New York, Oxford, 1978. 

Fakultat fur Mathematik, Universitat Regensburg, Universitatsstrasse 31, D-93053 
Regensburg, Germany 

E-mail address: helmut.abels@mathematik.uni-regensburg.de 

Institut fur Mathematik, Martin-Luther-Universitat Halle- Wittenberg, Theodor- 
Lieser-Str. 5, D-06120 Halle, Germany 

E-mail address: mathias.wilkeOmathematik.uni-halle.de 



